arXiv: astro-ph/0512271 v 1 11 Dec 2005 


On the Collapsar Model of Long Gamma-Ray Bursts: Constraints 

from Cosmic Metallicity Evolution 

N. Langer 1 2 3 and C.A. Norman 1,2,3 

ABSTRACT 

We explore the consequences of new observational and theoretical evidence 
that long gamma-ray bursts prefer low metallicity environments. Using recently 
derived mass-metallicity correlations and the mass function from SDSS studies, 
and adopting an average cosmic metallicity evolution from Kewley & Kobulnicky 
(2005) and Savaglio et al. (2005) we derive expressions for the the relative number 
of massive stars formed below a given fraction of solar metallicity, e, as function 
of redshift. We demonstrate that about 1 / 10th of all stars form with e < 0.1. 
Therefore, a picture where the majority of GRBs form with e < 0.1 is not in¬ 
consistent with an empirical global SN/GRB ratio of 1/1000. It implies that 
(1) GRB’s peak at a significantly higher redshift than supernovae; (2) massive 
star evolution at low metallicity may be qualitatively different and; (3) the larger 
the low-met allicity bias of GRBs the less likely binary evolution channels can be 
significant GRB producers. 

Subject headings: Gamma-ray bursts, Stars: Wolf-Rayet, Cosmic star formation, 
Galaxies 


1. Introduction 

It become clear in the last few years that long gamma-ray bursts are associated with 
the endpoints of massive star evolution. They occur in star forming regions at cosmological 
distances (Jakobsen et al. 2005), and are assiciated with supernova-type energies. The 
collapsar model explains gamma-ray burst formation via the collapse of a rapidly rotation 
massive iron core into a black hole (Woosley 1993). The short time scale of gamma-ray 
emission requires a compact stellar size, of the order of lightseconds. This constraint leaves 
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only massive Wolf-Rayet stars as possible progenitors. However, this poses a difficulty: 
Wolf-Rayet stars in the local universe are known to have strong stellar winds (Nugis et ah 
1998), which lead to a rapid spin-down (Langer 1998) — in agreement with the absence of 
signatures of rapid rotation in the Galactic Wolf-Rayet sample (Eenens 2004). 

It is the current understanding that the ratio of gamma-ray bursts to supernovae is 
about 1/1000, based on about 1-2 10 -6 observed bursts per supernova in the BATSE 
sample (Porciani & Madau 2001), and a beaming factor of ~ 500 (Frail et al. 2001; Yonetoku 
et al. 2005). This implies that about 1 out of 100 Wolf-Rayet stars produces a gamma-ray 
burst (van Putten 2004). These low values are thought to support the idea that rather exotic 
binary evolution channels might constitute the main evolutionary paths towards gamma-ray 
bursts (Podsiadlowski et al. 2004; Fryer and Heger 2005), corroborated by the two basic 
problems of single star models to produce collapsars (Petrovic et al. 2005): (1) the spin- 
down of the stellar core due to magnetic core-envelope coupling (Heger, Woosley & Spruit 
2005; Petrovic et al. 2005) — which is required to understand the slow rotation of young 
pulsars (Ott et al. 2005) and white dwarfs (Berger et al. 2005); and (2) the spin-down due 
to Wolf-Rayet winds mentioned above. 

However, recent single star models have overcome both problems for suitable initial 
conditions. The models by Woosley & Heger (2005) and Yoon & Langer (2005) avoid problem 
(1) by rapid rotation — which keeps the stars nearly chemically homogeneous and thus avoids 
the formation of a massive envelope — and problem (2) by choosing a low enough metallicity 
— which, according to recent evidence (Crowther et al. 2002; Vink & de Koter 2005), reduces 
the Wolf-Rayet mass loss rates. These are the first single star evolution models fulfilling the 
requirements of the collapsar model which are at the same time fully consistent with the 
slowly rotating stellar remnants in our Galaxy. However, they predict long GRBs only for 
metallicities of about Z/Z Q < 0.1. 

In this context, the growing empirical evidence that the long bursts indeed prefer a low- 
Z environment is remarkable. While indirect evidence from gamma-ray burst host galaxies 
is pointing towards low metallicity (Fynbo et al. 2003; Conselice et al. 2005; Fruchter et 
al. 2005), direct metallicity determinations yield sub-solar values down to l/100th of the 
solar metallicity (Gorosabcl et al. 2005; Chen et al. 2005; Starling et al. 2005). While 
the observational and the theoretical evidence for long gamma-ray bursts occurring at low 
metallicity needs further confirmation, we are motivated by the findings reported above to 
explore the consequences of such a possibility. 
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2. Method and results 

In order to quantify the amount of star formation at a given metallicity, we use the 
Schechter distribution function of galaxy masses and then substitute in it the well-defined 
mass-metallicity relation. We use the mass function and metallicity-mass correlation deter¬ 
mined in comprehensive SDSS studies and deep surveys. For the mass function there are 
fine studies by Cole et al. (2001), Bell et al. (2004) and Panter et al. (2004), and we use the 
latter here. For the mass-metallicity relation we use the classic studies of Tremonti et al. 
(2005) and Savaglio et al. (2005). From Equation (1) of Panter et al. (2004) 


*(M) = #, (41) V" 7 "-’, (1) 

where a = —1.16 and <f>* = 7.8 x 10~ 3 h 3 Mpc~ 3 . The fraction T(M) of the mass density in 
galaxies with a mass less than M is then 


T(M) 


f 0 M M$(M)dM 
f 0 °° M$(M)dM 

f (a + 2, (M/M*)) 

r(a + 2) 


( 2 ) 


where f and T are the incomplete and complete gamma function. We then choose the galaxy 
mass-metallicity relation of the form: M/M* = K^Z/Zq) 13 where K and /3 are constants. 
For simplicity we choose to use the linear bisector fit to the mass metallicity relation derived 
by Savaglio et al. (2005) at redshift, z = 0.7. This is parallel to the quadratic fit of Tremonti 
et al. (2005) at low masses and metallicities that we consider here. The form we derive is 
very close to /3 — 2 and K = 1 and we assume these values. This gives a pleasing Gaussian 
form to the metallicity function which may be of more general applicability. We now derive 
the fractional mass density belonging to metallicities below metallicity Z at a given redshift 
z as 


^ (Z_\ = f(q + 2,(Z/Z Q flO°- 15 ^) 
\Zq) T(a + 2) 


( 3 ) 


where we used the average cosmic metallicity scaling as = —0.15 dex per unit redshift 
from Kewley & Kobulnicky (2005, 2006). We used this scaling since the mean metallicity is 
given by: 


(Z) 
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and thus the mean metallicity is linearly proportional to Z Q , (Z) = constant x Z Q . The 
ratio of the gamma-functions in the above equation gives 0.80 for (3 — 5, and 0.73 for (3 = 2, 
assuming a = —1.16 in both cases. Our scaling simply gives self-consistently (Z) = Z e 10~ 7Z 
where from Kewley and Kobulnicky(2005, 2006) we have 7 = 0.15 that we adopt here. The 
7 value found by Savaglio et al. (2005) is approximately 7 ~ 0.3 — 0.4 (S.Savaglio, private 
communication). 

This ansatz includes various simplifications. Because the parameter we need here is a 
ratio of two moments of the mass function the evolution of the normalisation of the mass 
function (the cosmic stellar mass density) cancels. The remaining evolutionary term is 
essentially the mean cosmic metallicity evolution. Since we are dealing here with massive 
star formation from the interstellar medium and since the metallicity constraint is from 
the physics of line-driven winds we use the scaling with redshift derived from emission-line 
studies by Kewley and Kobulnicky (2005, 2006) for [O/H]. The study of Gallazzi et al. (2005) 
derives stellar metallicities and gives in their Figure 6 the distribution of stellar metallicity 
in galaxies at low redshift. As they note, in general, the stellar metallicity is always lower 
than the gas-phase metallicity and our use of Kewley and Kobulnicky’s gas-phase evolution 
is conservative. We are assuming in these star-forming galaxies that the evolving metallicity 
distribution we have derived also applies to the ISM in which the stars are forming. 


For our fiducial numbers a = —1.16, (3 = 2 and K = lwe thus obtain 





f (0.84, (Z/Z o ) 2 10°- 30 ^) 
r(0.84) 


where T(0.84) ~ 1.122. 


( 5 ) 


The result of folding T with the total star formation rate history is shown in Figure 1. 
Here, the total star formation rate rgFR^) is derived from a 4th-order polynomial fit to the 
data presented in Bouwens et al. (2004). 


In order to be able to compare with observations, we convolve the fraction T of stars 
born with a metallicity of less than Z Q e with the comoving volume element of the Friedman- 
Robertson-Walker metric. The quantity d m (z ) used below is the proper motion distance 
for our cosmological parameters. We are using the standard cosmological parameters of 
H 0 = 70 krns -1 Mpc -1 , Hm = 0.3, and Ha = 0.7. 


Including the time dilation effect due to redshift that slows the rates with redshift, the 
unbiased observed rate of core collapse supernovae from stars with a metallicity below Z Q e 
then becomes 


Rsv(e) — /sn4vt 
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where 

plOO Mq 

/sn = / atp(rn)drn ~ 0.0074 (7) 

J 8 Mq 

is the number of core collapse events per solar mass of stars formed, where a = 1 M 0 / J ip(m)dm, 
(p(m) is the Salpeter initial stellar mass function, and /sn is the number of core collapse events 
per solar mass of star formation. Figure 2 shows the resulting supernova rates, and their 
integrated value, for various values of e. 

In Table 1, we derive the preceived global and local ratios of low-metallicity to all su¬ 
pernovae, for various metallicity thresholds and for two galaxy mass-metallicity exponents. 
From these, we compute formation rates of low-metallicity black holes by assuming super- 
novae to form from stars above 8M 0 , but black holes to form only from stars above 30 M 0 , 
which — for a Salpeter IMF — gives a BH/SN ratio of 0.14. Since black hole formation is 
only one out of three criteria for GRB production within the collapsar frame, our derived 
numbers give an upper limit to the GRB production rate from low metallicities, independent 
of the stellar evolution scenario considered. For £ = 0.1, this corresponds to less than 22 
GRBs per 1000 SNe globally in the universe, and to less than 3/1000 at low redshift. This 
is what stars with a metallicity of Z Q /10 and below can provide maximally. Table 1 shows 
further that the constraint of achieving 1 GRB per 1000 SNe globally in the universe results 
in about 1 GRB per 10 000 SNe in the local universe. 

We note that these numbers seem to be remarkably insensitive to details of the star 
formation history. The last row in Table 1 is computed by using the star formation history 
shown by Firmani et al. (2005), which is similar to the one used here up to redshift 2, but 
which remains at the top level until about z = 6 (cf. their Fig. 2). However, Table 1 also 
shows that the redshift at which the GRB rate peaks does depend on the star formation 
history and can thus not yet be predicted reliably. 


3. Discussion and conclusions 

To produce a GRB, Wolf-Rayet stars at core collapse are required to have sufficient 
angular momentum. Stellar evolution models which include magnetic fields predict too 
slowly-rotating cores for models which develop an extended, massive envelope after the main 
sequence. Current evolutionary models that include rotation predict extended envelopes for 
the vast majority of massive stars in the Galaxy or Magellanic Clouds — in agreement with 
the number of blue and red supergiants (Maeder & Meynet 2001, 2005). Only the fastest 
rotators are thought to be able to avoid extended envelopes, for metallicities below about 
Z Q /10, or £ = 0.1 (Yoon & Langer 2005; Woosley & Heger 2005). 
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Most importantly, the numbers derived for £ = 0.1 in Table 1 appear not to be in 
conflict with observations. Per 1000 supernovae in the universe, 160 are predicted to occur 
from stars with Z > Z 0 /1O, out of which 22 would produce a black hole. Thus, producing 
one GRB per 1000 supernovae globally in the universe (cf. Sect. 1) seems possible. 

On the other hand, a large fraction of the 22 black holes may be born without producing 
a GRB: not all of them may occur in WR stars but rather in more extended stars (Maeder 
& Meynet 2005), and the most massive ones would lose too much angular momentum in 
a wind, even for metallicities as low as Z > Z Q j 10 (Yoon & Langer 2005). While within 
the chemically-homogeneous-evolution scenario for GRB formation (Yoon & Langer 2005; 
Woosley & Heger 2005) a GRB/BH fraction of 1/20 can certainly be obtained (Yoon & 
Langer 2006), this is therefore unlikely for exotic binary channels for GRB production — i.e. 
for channels through which only a small fraction of stars of any initial mass evolves. Clearly, 
the more the long GRBs are confined to low metallicities, the more unlikely it is that binary 
evolution is needed to explain the majority of events. 

The empiric cosmic GRB to SN ratio of about 1/1000 (cf. Sect. 1) can not directly rule 
out more extreme values of e, i.e. £ = 0.1 — 0.01 (cf. Table 1); in fact, for e = 0.01, about 
10% of all massive stars with Z < Z Q j 100 would need to produce a GRB. However, it would 
imply that the formation of every black hole would be accompanied by a GRB. Furthermore, 
the GRB rate would peak only at a redshift of about z — 10. A value of e = 0.01 appears 
thus unlikely. Furthermore, our models with j5 = 5 produce such a small local GRB/SN 
ratio that they seem to be ruled out. 

We find that a restriction of GRBs to low metallicities (Z < Z Q / 10, i.e. £ = 0.1) has 
the following consequences: 

• GRBs do not follow star formation in an unbiased manner. For example, for an overall 
star formation rate which predicts a preceived SN peak at a redshift of zsm — 1.8 we 
find, for £ = 0.1, that the GRB rate peaks at a redshift of zgrb — 3.2 (Fig. 2 and 
Tab. 1; see also Firmani et al. (2005) 

• Local massive galaxies, like our Milkyway, are not expected to host long GRBs. The 
last long GRB in our Galaxy should have occurred several gigayears ago. 

• The global and local GRB to SN ratios appear to be insensitive to the details of 
the cosmic star formation history, while the redshift of the peak GRB rate can vary 
appreciably (Tab. 1). 

• The local GRB/core-collapse ratio is much smaller than the one obtained from aver¬ 
aging over the universe; i.e., by one order of magnitude for e = 0.1 (Tab. 1). 
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• We obtain the expected result that the number of massive stars in the universe with 
a metallicity below a critical value £ does roughly scale with e. I.e., For £ = 0.1, we 
find a ratio of low-metallicity (Z < Z Q /10) to total global supernova rates of 0.16. 
Further, the probability that a randomly chosen burst has a metallicity of Z Q /100 is 
about 10%. 

• We derive the most likely redshift for gamma-ray bursts of specified metallicity for unbi¬ 
ased observations. For example, for a metallicity of Zq/ 100 - as found for GRB050730 

, we find GRBs to occur most likely at redshifts of z > 6. Locally, the ratio of GRBs 
with a metallicity of Z Q /100 to all GRBs is about 0.02 (Tab. 1). 

• The larger the low-metallicity bias of long GRBs, the less likely can binary scenarios 
explain the major fraction of them. 

• A confirmation of the low-metallicity bias of long GRBs to values of the order of £ = 0.1 
would imply that fast rotation may be much more common at low metallicity among 
massive stars. 


We are grateful to Bram Achterberg, Andrew Fruchter, Karl Glazebrook, Max Petini, 
Sandra Savaglio, Rosemary Wyse and Stan Woosley for very interesting and useful discus¬ 
sions. C.N. is very pleased to thank the Astronomical Institute in Utrecht for stimulating 
hospitality. 
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Fig. 1.— Formation rate of stars with a metallicity below Z Q /10, (e = 0.1) for a fiducial 
galaxy mass-metallicity relation M ~ Z 13 for f3 = 2 (blue curve) and 0 — 5 (red curve). 
These curves have been obtained by multiplication of the fit to the total star formation 
rate as function of redshift (weak green line) obtained from the empirical data presented 
by Bouwens et al. (2004; points with error bars) with the fraction T of stars born with a 
metallicity below a specified value. 







fobs SN, low-Z SN rate dz obs. SN, low-Z SN rates per dz 
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Fig. 2.— Normalised differential and cumulative core collapse supernovae rates as function 
of redshift as preceived by an unbiased observer (cf. Eq. 6), for all stars, and for stars with 
metallicities less than Z 0 /3 and Zq/ 10. 
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Table 1: For various upper metallicity limits e, and galaxy mass-metallicity exponents (3\ 
the number of supernovae with Z/Z & < e (or “low -Z supernovae”) per 1000 supernovae in 
the universe as seen from an unbiased observer (cf. Eq. 6, and Fig. 2), the corresponding 
number of “low -Z black holes” per 1000 supernovae, both numbers at redshift z — 0, the 
local (z = 0) number of GRBs per 1000 supernovae required to have a global GRB/SN ratio 
of 1/1000, and the redshifts at which the SN and GRB-rates are preceived as maximum by 
an unbiased observer (cf. Fig. 2). For the last row, a different underlying star formation 
history has been assumed (see text). 
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